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We utilize a superconducting stripline resonator containing a dc-SQUID as a strong intermod- 
ulation amplifier exhibiting a signal gain of 25 dB and a phase modulation of 30 dB. Studying 
the system response in the time domain near the intermodulation amplification threshold reveals a 
unique noise-induced spikes behavior. We account for this response qualitatively via solving numer- 
ically the equations of motion for the integrated system. Furthermore, employing this device as a 
parametric amplifier yields a gain of 38 dB in the generated side-band signal. 

PACS numbers: 05.45.-a, 85.25.Dq, 84.40.Dc 
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The field of solid-state qubits and quantum informa- 
tion processing has received considerable attention dur- 
ing the past decade P, 0, H, |3| and has successfully 
demonstrated several milestone results to date [1, H, 0, H, 
B [13, [HI- However, one of the fundamental challenges 
hindering this emerging field is noise interference, which 
either screens the output signal or leads to quantum state 
decoherence [1, H, [12, IH, [13, [IB]- Hence, this may ex- 
plain to some extent the renewed interest exhibited re- 
cently by the quantum measurement communit y in the 
field of phase-sensitive amplifiers [H, [13, [H, [H 113, [H • 
This interest is driven essentially by two important prop- 
erties of these amplifiers: (1) Their capability to amplify 
very weak coherent signals; (2) Their ability to squeeze 
noise below the equilibrium level by means of employ- 
ing a homodyne setup and phase control. These prop- 
erties are expected to be highly beneficial to the area 
of quantum communication |22l . |23|| and to t he g enera- 
tion of quantum squeezed states [ll, [13, [H, [IBl- Ad- 
ditional interest in these high-gain parametric amplifiers 
arises from the large body of theoretical work predict- 
ing photon-generation from vacuum via the dynamical 
Casimir effect [26', '27?], which can be achieved by em- 
ploying an appropriate parametric excitation mechanism 

In this present work we study a superconducting 
stripline microwave resonator integrated with a dc- 
SQUID [13, HH, HI]. The paper is mainly devoted to 
a novel amplification mechanism in which the relatively 
large nonlinear inductance of the dc-SQUID is exploited 
to achieve large gain in an intermodulation (IM) con- 
figuration. We provide theoretical evidence to support 
our hypothesis that the underlying mechanism responsi- 
ble for the large observed gain is metastability of the dc- 
SQUID. In addition, at the end of the paper we demon- 
strate another amplification mechanism, in which the de- 



* Electronic address: ' baleegh@tx. technion. acTiT] 



pendence of the dc-SQUID inductance on external mag- 
netic field is exploited to achieve parametric gain j20j . 

The device (see Fig. [T] (a)) is implemented on a high- 
resistivity 34 mm X 30 mm X 0.5 mm Silicon wafer coated 
with a 100 nm thick layer of SiN. As a preliminary step, 
thick gold pads (300 nm) are realized at the periphery of 
the wafer. Subsequently, e-beam lithography is applied 
in which both the macroscopic resonator and the micro- 
scopic Josephson gaps are written. Following a develop- 
ing stage, a two angle shadow evaporation of Aluminium 
[33| -with an intermediate stage of oxidation- implements 
the resonator as well as the two Al/AlO^^/Al Josephson 
junctions comprising the dc-SQUID. The total thickness 
of the Aluminium evaporated is 80 nm (40 nm at each an- 
gle). Finally, a lift-off process concludes the fabrication 
of the integrated system. 

In general, IM generation is a nonlinear phenomenon 
which is associated often with the occurrence of nonlin- 
ear effects in the resonance curves of a superconduct- 
ing resonator [H, [H, [H, [13, [111 . In Fig. [1 we show a 
transmission measurement of the resonator response, ob- 
tained using a vector network analyzer at the resonance 
while sweeping the frequency upwards. The measured 
resonance resides at /o = 8.219 GHz which corresponds 
to the third resonance mode of the resonator having an 
anti-node of the rf-current waveform at the dc-SQUID 
position. Similar nonlinear effects in the transmission re- 
sponse have been measured at the first mode (~ 2.766 
GHz) as well. 

As can be seen in the figure, at excitation powers Pp be- 
low Pc = —60.3 dBm (at which nonlinear effects emerge) 
the resonance curve is linear and Lorentzian. As the in- 
put power is increased abrupt jumps appear at both sides 
of the resonance (see black line in Fig. [2]) accompanied by 
frequency hysteresis loops. In addition, as one continues 
to increase the input power the resonance curves become 
shallower, broader and less symmetrical. Such power de- 
pendency can be attributed to the nonlinearity of the 
dc-SQUID inside the resonator which increases consid- 
erably as the amplitude of the driving current becomes 
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FIG. 1: (Color) The device and IM setup, (al) A photograph 
of the device consisting of a circular resonator with a circum- 
ference of 36.4 mm and a line width of 150 ^m, a dc-SQUID 
integrated into the resonator and a flux-line employed for driv- 
ing rf-power and flux-biasing the dc-SQUID. (a2) An electron 
micrograph displaying the dc-SQUID incorporated in the res- 
onator and the adjacent flux-line. The area of the dc-SQUID 
fabricated is 39 fim x 39 /xm, while the area of each junction 
is about 0.25(/im)^. The self-inductance of the dc-SQUID 
loop is Ls = 1.3 ■ 10"^'' H and its total critical Josephson cur- 
rent is about 2/c = 28 fj, A. (a3) A micrograph of one of the 
dc-SQUID junctions, (b) A basic IM setup. The pump and 
the signal designated by their angular frequency Up and tOs 
respectively are combined by a power combiner and fed to 
the resonator. The field at the output is amplified using two 
amplification stages and measured using a spectrum analyzer 
(path 1). A dc current was applied to the flux- line in order to 
flux-bias the dc-SQUID. Path (2) at the output corresponds 
to a homodyne setup employed in measuring IM. 



comparable to the critical current. 

The basic IM scheme used is schematically depicted in 
Fig. [1] (b). The input field of the resonator is composed 
of two sinusoidal fields generated by external microwave 
synthesizers and superimposed using a power combiner. 
The applied signals have unequal amplitudes. One, re- 
ferred to as the pump, is an intense sinusoidal field with 
frequency /p, whereas the other, referred to as the sig- 
nal, is a small amplitude sinusoidal field with frequency 
fp+S, where S represents the frequency offset between the 
two signals. Due to the presence of a nonlinear element 
such as the dc-SQUID integrated into the resonator, fre- 
quency mixing between the pump and the signal yields 
an output idler field at frequency fp — 5. Thus the output 
field of the resonator, measured by a spectrum analyzer, 
consists mainly of three spectral components, the trans- 
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FIG. 2: (Color) Transmission response curves of the resonator 
at its third resonance (8.219 GHz) as a function of input power 
and frequency. The resonance curves exhibit nonlinear effects 
at Pp > Pc- The loaded quality factor for this resonance in the 
linear regime is 350. The black line maps the jumps appearing 
at both sides of the resonance curves. The inset exhibits the 
dependence of the transmitted pump signal on the applied 
magnetic flux threading the dc-SQUID loop during an IM 
measurement. 



mitted pump, the transmitted signal and the generated 
idler. The IM amplification in the signal, idler and the 
noise is obtained, as shown herein, by driving the dc- 
SQUID to its onset of instability via tuning the pump 
power. 

In Fig. [3] we show a typical IM measurement result. In 
this measurement the pump was tuned to the resonance 
frequency fp — 8.219 GHz. The signal power Ps was set 
to —100 dBm, and its frequency offset 6 to 5 kHz. As the 
pump power injected into the resonator is increased, the 
nonlinearity of the dc-SQUID increases and consequently 
the frequency mixing between the pump and the signal 
increases as well. At about Pc the dc-SQUID reaches 
a critical point at which the idler, the signal and also 
the noise floor level (within the frequency bandwidth) 
undergo a large simultaneous amplification. The idler 
gain measured relative to the injected signal power at 
the resonator input is about 13 dB (see inset of Fig. [3]). 
The same result is obtained as well, as one calculates the 
power ratio of the generated signal at the output port 
of the resonator to the input signal while accounting for 
the losses and amplifications of the elements along each 
direction. 

In order to show that this IM amplifier is also phase 
sensitive we applied a standard homodyne detection 
scheme as schematically depicted in Fig. [T] (see path 
(2)), in which the phase difference between the pump and 
the local oscillator (LO) at the output -having the same 
frequency- can be varied. In such a scheme the pump 
is down-converted to dc, whereas both signal and idler 
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FIG. 3: (Color) A spectrum power measured by a spectrum 
analyzer during IM operation as a function of increasing pump 
power while applying external flux $ ~ $o/4. The spectrum 
taken at a constant frequency span around fp was shifted 
to dc for clarity. At the vicinity of Pc the system exhibits 
a large amplification. A cross-section of the measurement 
taken along the pump power axis is shown in the inset. The 
red line (light grey) and the blue line (dark grey) exhibit the 
corresponding gain factors of the transmitted signal and the 
idler respectively. 



are down-converted to the same frequency S. The largest 
amplification is obtained when the LO phase {(f>Lo) is ad- 
justed such that the signal and idler tones constructively 
interfere. Shifting (j)Lo by tt/2 from the point of largest 
amplification results in destructive interference, which in 
turn leads to the largest de-amplification. 

The IM measurement results, obtained using the ho- 
modyne setup while flux-biasing the dc-SQUID with 
about half flux-quantum, arc shown in Figs. |3]and[ni Fig- 
ured] exhibits the signal gain (blue line) versus increasing 
applied pump power. As can be seen from the figure the 
signal gain assumes a peak of 25 dB for Pp :^ Pc- More- 
over, the response of the resonator at the pump frequency 
(green line) -measured simultaneously using a voltage 
meter connected in parallel to the spectrum analyzer- 
is drawn as well for comparison. The sharp drop in the 
region Pp ~ Pc, coinciding with the amplification of the 
signal, indicates that the transmitted pump is depleted 
and power is transferred to other frequencies. 

Fig. ^exhibits a periodic dependence of the signal gain 
on 0iOj having a gain modulation of about 30 dB peak 
to peak between the ampliflcation and de-amplification 
quadratures. It is worthwhile mentioning that in this 
measurement not only the signal shows a periodic de- 
pendence on <f>Lo but also the noise floor which exhibits 
up to 20 dB modulation (measured at fp + l.b5). 

Another interesting aspect of this device is revealed by 
examining its response in the time domain while applying 
an IM measurement using the homodyne setup shown in 
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FIG. 4: (Color) A signal gain (blue line) measured as a func- 
tion of increasing pump power. The measurement was taken 
using a homodyne scheme. The signal displays a gain of about 
25 dB at Pc. The corresponding dc component of the homo- 
dyne detector output (green line) was measured simultane- 
ously using a voltage meter in parallel to the spectrum ana- 
lyzer. In this measurement Ps = —110 dBm and 6 = 5 kHz. 
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FIG. 5; A periodic dependence of the signal gain on the LO 
phase difference at the vicinity of Pc. The signal exhibits a 
large amplification and de-amplification at integer and half 
integer multiples of tt respectively. The phase modulation 
dependency shows up to 30 dB peak to peak amplitude. 



Fig. [T] (see path 2). This is achieved by connecting a 
fast oscilloscope in parallel to the spectrum analyzer at 
the output. A few representative snapshots of the device 
behavior measured in the time domain are shown in Fig. 
El where each subplot corresponds to a different applied 
pump power. In subplot (a) the pump power is about 
3 dBm lower than the critical value and it displays a 
constant voltage level. As the pump power is increased 
(subplot (b)) separated spikes start to emerge. Increasing 
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FIG. 6: Snapshots of the resonator response measured in the 
time domain during an IM operation, (a)-(d) display the 
resonator response at pump powers lower than the critical 
value, (e) At the critical value, (f) Exceeding the critical 
value. 



the power further causes the spikes to start bunching with 
one another and as a consequence forming larger groups 
(subplots (c) and (d)). This bunching process reaches an 
optimal point at the pump power corresponding to the 
peak in the IM gain (subplot (e)). Above that value, such 
as the case in subplot (f), the spikes become saturated 
and the gain drops. 

Such behavior in the time domain reveals also the un- 
derlying mechanism responsible for the amplification. As 
it is clear from the measurement traces the rate of spikes 
strongly depends on Pp at the vicinity of Pc, thus adding 
a small signal at fp + S (which is effectively equivalent to 
applying an amplitude modulation of the pump at fre- 
quency S), results in modulation of the rate of spikes at 
the same frequency, and consequently yields a large re- 
sponse at the signal and the idler tones. 

In an attempt to account for the unique temporal re- 
sponse presented in Fig. [6l we refer to the two equations 
of motion governing the dc-SQUID system given by [3§| 
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where $0 is the flux quantum, Rj is the shunt resis- 
tance of the junctions, Ic is the critical current of each 
junction, iJi, 62 are the gauge invariant phase differences 
across junctions 1 and 2 respectively, Ini,In2 are equi- 
librium noise currents generated in the shunt resistors 
having -in the limit of high temperature- a spectral den- 
sity of Sj^ = AksT/Rj, Pc is a dimensionless parameter 



defined as f3c = 27r/ci?jC,//$o, where Cj is the junction 
capacitance and U is the potential energy of the system 
which reads 



U = -^-{d^ + S2) + ^ 

— Ic (cos 61 + COS S2 ) , 
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where $ is the applied magnetic flux, / is the bias current 
flowing through the dc-SQUID, and /3l is a dimensionless 
parameter defined as (3^ = 2LsIc/^o- While the circu- 
lating current flowing in the dc-SQUID loop is given by 
/circ (<) - Ic [Si - 52 - 27r$/$o) /tt/Jl. 

Furthermore, in order to account for the resonator re- 
sponse as well, we make two simplifying assumptions: 
(1) We model our resonator as a series RLC tank os- 
cillator characterized by an angular frequency luq = 
1 / VLC = 27r • 8.219 GHz, and a characteristic impedance 
Zq ~ -y/ L/C, where L and C are the inductance and the 
capacitance of the resonant circuit respectively; (2) We 
neglect the mutual inductance that may exist between 
the inductor and the dc-SQUID. 

Under these simplifying assumptions one can write the 
following equation of motion for the charge on the capac- 
itor denoted by q {t) 

— +Rq + UJoZoq + Vsq + Vin + U„ = 0, (3) 

UJQ 

where is a voltage noise term having -in the limit of 
high temperature- a spectral density of Sy^ = '^RksT, 
Vi-n (<) = Re (Vbe~'"p*) is an externally applied sinu- 
soidal voltage oscillating at the pump angular frequency 
LOp and having a constant voltage amplitude Vq. Whereas, 

^sq = $0 i^i + 52^ 1^1^ designates the voltage across the 
dc SQUID. Using these notations the bias current in Eq. 
dU) reads I = q. 

Finally in order to relate the observed output signal in 
Fig. [H] at the output of the homodyne scheme to the 
fast oscillating solution g, we f irst ex press the charge 
on the capacitor as q(t) = \/2h/Zo Re [A (t) e"*""*] , 
where A (t) is a slow envelope function, second, we em- 
ploy the following input-output relation [13] Vout (t) = 
Vq — iy/32ZohjfA (t), where 71 is the coupling constant 
between the resonator and the feedline, in order to ob- 
tain the field at the output port, third, we account for 
the phase shift by evaluating the following expression 
2Vlo Re [Vl^ut (0 6"^^°] , where Vlo corresponds to the 
amplitude of the LO. 

By integrating these stochastic coupled equations of 
motion numerically, while employing device parameters 
which are relevant for our case, one finds that the ob- 
served temporal behavior of the system can be qualita- 
tively explained in terms of noise-induced jumps between 
different potential wells forming the potential landscape 
of the dc-SQUID which is given by Eq. As a con- 

sequence of applying the voltage Vin to the integrated 
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FIG. 7: A simulation result. The output field exhibits spikes 
in the time domain similar to the spikes shown in Fig. (6] The 
spikes occur whenever the system locus jumps from one well 
to another due to the presence of stochastic noise while driv- 
ing the system near a critical point. In this simulation run 
the system was simulated over a period of 1000 time cycles of 
pump oscillations (one of the constraints set on the simulation 
was maintaining a reasonable computation time). The param- 
eters that were employed in the simulation (same as experi- 
ment) are ujq = ujp ^ 2tt ■ 8.219 GHz, /3l = 3.99, Rj = 9Afl, 
4" = 'E'o/S and (I>lo = ti". The rest of the parameters were set 
in order to reproduce the measurement result: 71 = 10^ Hz, 
R = 3.3 f2, (3c = 5.86 and ^0 = 10 f2 (corresponding experi- 
mental values are: 71 ~ 2.4 • 10^ Hz, Zo ~ son, Cj was not 
measured directly, a capacitance on the order of 0.7 pF was 
assumed) . 



system, the potential landscape of the dc-SQUID os- 
cillates at the pump frequency, and its oscillation am- 
plitude grows with the amplitude of the incoming volt- 
age. However, inter- well transitions of the system become 
most dominant as the oscillation voltage reaches a criti- 
cal value at which the potential landscape of dc-SQUID 
becomes tilted enough -due to the current flowing in the 
system- in order to allow frequent noise-assisted escape 
events from one well to another or across several wells, 
which in turn cause a voltage drop to develop across the 
dc-SQUID and induce jumps in the circulating current. 
Such hopping of the system state is manifested as well 
in the homodyne output field response shown in Fig. [3 
which in a qualitative manner mimics successfully the 
main temporal features shown in Fig. [Sj 

It is also clear from this inter-well transitions model 
that the ampliflcation in this case is different from the 
so-called Josephson bifurcation amplifier [l^, as in the 
latter case the system is confined to only one well and the 
bifurcation occurs between two oscillation states having 
different amplitudes. 

Moreover, just as in the recent experiment by Ya- 
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FIG. 8: (Color) (a) A homodyne setup employed in measur- 
ing parametric excitation, (b) Large gain measured in a para- 
metric excitation experiment corresponding to the generated 
three harmonics (at S — 5 kHz, 2S, 3S respectively) and the 
noise spectral density at 1.55. The gain was measured as a 
function of increasing pump power applied at 2fo+5. The sig- 
nal power applied to the resonator at /o was Ps = —80 dBm. 
In this measurement a maximum gain of 38 dB is achieved at 
the first harmonic. 



mamoto et al. j20| we have additionally employed our 
device as a parametric amplifier. To this end, we have 
used the parametric excitation scheme exhibited in Fig. 
[5] (a) in which the pump and signal tones are applied to 
different ports. The main rf signal (pump) is applied to 
the flux-line at the vicinity of twice the resonance fre- 
quency of the resonator 2/o -I- <5 (i5 = A/27r = 5 kHz) and 
2/0 — 16.438 GHz does not coincide with any resonance 
of the device. Whereas, the signal, being several orders of 
magnitude lower than the pump, is fed to the resonator 
port at /o and its main purpose is to probe the system 
response. 

In Fig. [5] (b) we exhibit a parametric excitation mea- 
surement result obtained using this device at $ ~ 0.6$o- 
In this result there is evidence of one of the charac- 
teristic fingerprints of a parametric amplifier: the exis- 
tence of an excitation threshold above which there is a 
noise rise and an abrupt amplification of the harmonic at 
/o + (5 (which results from a nonlinear frequency mixing 
at the dc-SQUID). As can be seen in the figure at about 
—40 dBm the first harmonic generated at 5 and the two 
higher-order harmonics at 25, 3(5 get amplified consider- 
ably up to a maximum gain of 38 dB measured at the 
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first harmonic. Also in a separate measurement resiilt 
(not shown here) the first harmonic has been found to 
display ~ 20 dB peak to peak modulation as a function 
of external magnetic field. 

In conclusion, in this work we have designed and fab- 
ricated a superconducting striplinc resonator containing 
a dc-SQUID. We have shown that this integrated system 
can serve as a phase sensitive amplifier. We have studied 
the device using IM measurement and parametric exci- 
tation. In both schemes the device exhibited distinct 
threshold behavior, strong noise rise and large amplifica- 
tion of coherent side-band signals generated due to the 
nonlinearity of the dc-SQUID. In addition, we have inves- 
tigated the system response in the time domain during 
IM measurements. We have found that in the vicinity of 
the critical input power the system becomes metastable 
and consequently exhibits noise-activated spikes in the 
transmitted power. We have shown that this kind of 



behavior can be explained in terms of noise-assisted hop- 
ping of the system state between different potential wells. 
We have also demonstrated that the main features ob- 
served in the time domain can be qualitatively captured 
by solving the equations of motion for the dc-SQUID in 
the presence of rf-current bias and stochastic noise. Such 
a device may be exploited under suitable conditions in a 
variety of intriguing applications ranging from generat- 
ing quantum squeezed states to parametric excitation of 
zero-point fluctuations of the vacuum. 
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